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We analyze the outcomes of an extended-Drude-model approach to the optical spectra of pnic- 
tides, where the multiband nature of the electronic excitations requires a careful analysis of the 
role of interband processes in the optical conductivity. Through a direct comparison between model 
calculations of the intraband optical spectra and experimental data, we show that interband tran- 
sitions, whose relevance is shown by first-principle calculations, give a non negligible contribution 
already in the infrared region. This leads to a substantial failure of the extended-Drude-model 
analysis on the measured optical data without subtraction of interband contributions. 



Optical studies are a useful experimental probe to an- 
alyze the interactions at play in different classes of cor- 
related materials, 1 for which it would be desirable to 
devise a common theoretical scheme. A typical exam- 
ple is the so-called extended-Drude-model (EDM) anal- 
ysis, where the experimental optical conductivity is an- 
alyzed in terms of a Drude-like model with a frequency- 
dependent inverse lifetime and effective massi This ap- 
proach is well justified for a single-band system inter- 
acting with a continuum of bosonic excitations (such for 
instance phonons in a metal) where it can provide useful 
information on the relevant collective modes interacting 
with the electronic particle-hole excitations, 2 The EDM 
has been widely used also for strongly correlated materi- 
als, such as the cuprate superconductors; here, intraband 
excitations are well separated from the optical interband 
transitions, which pose a natural cut-off (typically of or- 
der of 4000 cm" 1 ) to the applicability of the EDM itself 
in these compounds 

After the discovery of superconductivity in iron-based 
superconductors an intense experimental and theoretical 
research has been devoted to the EDM analysis in these 
materials as welli^— A general outcome of this analysis 
is a relatively large frequency-dependent in-plane inverse 
lifetime t _1 (oj) ~ lu,^^ that, within the EDM, is inter- 
preted in terms of a strong-coupling regime (A ~ 3 — 4) 
for the relevant bosonic excitations located around 20- 
60 meV. These are usually identified with the spin fluc- 
tuations between the hole and electron pockets of the 
Fermi surface, that from the very beginning have been 
suggested as the most promising candidates for the su- 
perconducting pairing. 9 

There are however two main open issues concerning 
these results. On one hand, a strong coupling to these 
low-energy bosonic modes would imply also an effective 
mass at low energy much larger than that measured by 
other probes such as specific-heafeii or de Haas- van 
Alphen (dHvA)ji^ On the other hand, the use of the 
EDM implicitly assumes that the interband transitions 
in pnictides become relevant only above a threshold fre- 



quency of about ~ 2000 cm -1 , even though some authors 
suggested that they are present at lower energies J£r— 
Low-energy interband transitions may invalidate the con- 
clusions of the EDM analysis, introducing spurious ef- 
fects. 

In this paper, we quantify the relative role of inter- 
band and intraband transitions in the EDM analysis of 
the in-plane optical conductivity in pnictides, focusing 
on LaFePO where no magnetic transition occurs and 
where the coupling to spin fluctuations within a multi- 
band Eliashberg approach has been already discussed in 
the context of different physical properties! 16 ! 17 Through 
a direct comparison with the experimental optical data 
of Ref. [4] we show that, for the physical range of the 
coupling of electrons with spin-fluctuations, the intra- 
band contribution alone is not sufficient to account for 
the large inverse lifetime obtained by the EDM analysis. 
In this context we also show that low-energy interband 
transitions, individuated through first-principles calcu- 
lations, play a dominant role here already at very low 
energy (ui > 500 cm -1 ), invalidating the conclusion of a 
simple EDM analysis of the data. Such a result, along 
with the similarity of the band structures and optical 
spectra of Fe pnictides, suggests that a new approach is 
needed in order to analyze the optical data of pnictides. 

The EDM analysis is based on the fact that many- 
body interactions modify the Drude formula for the op- 
tical conductivity by introducing a frequency dependence 
into the inverse lifetime t -1 (oj) and a mass enhance- 
ment factor m*(ui)/mb, where mb is the band mass of 
the carriers^ Deconvoluting the real <7i(u>) and imagi- 
nary 02 (w) part of the complex optical conductivity, one 
can extract: 

H = 4^ 2 M+a 2 2 M' (1) 
m*(u}) = uj 2 p a 2 (uj) 
m A-klo (j\{lu) + o\{uS) ' 

where up is the plasma frequency, that is usually esti- 
mated integrating the real part of the conductivity up to 
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FIG. 1: (color online) Comparison between experimental data 
(filled circles, after Ref. |j) and theoretical calculations of 
the intraband contribution (lines) for the optical conductiv- 
ity (panel a) and the corresponding t _1 (cj) as extracted by 
an EDM analysis (panel b). Different values of the average 
coupling strength are here considered (see text). 



a cut-off Lu c - The underlying idea is that the data below 
lu c represent intraband optical transitions of interacting 
charge carriers (with spectral weight proportional to uj 2 p ), 
while the spectrum above ui c represents interband tran- 
sitions. In the case of pnictides ui c is usually assumed 
around 2000 - 3000 cm^ii 

Experimental data for o~i(uj) and for the resulting 
r _1 (w) in LaFePO, from Ref. 0, are shown in Fig. [1] 
as blue symbols. In particular, the analysis of r _1 (w) 
clearly reveals two distinct regimes which are also visible 
in tri(w): a low-energy one below ~ 500 cm" 1 , where 
t _1 (cj) is roughly constant (t -1 (w) rs 400cm ) and 
which corresponds to a very narrow Drude-like contri- 
bution in cti(cj); and a higher-energy regime [uj > 500 
cm -1 ) where oi(lo) is almost constant and T -1 (k;) has 
a steady increase, with t" 1 (oj) ~ ui. Note that scatter- 
ing by impurities leads to a constant t _1 (cl>), while scat- 
tering by retarded interactions is usually reflected in an 
increase of t~ 1 (cj) up a saturation regime at frequencies 
larger than the typical energy scale of the interaction. 2 
Since phonon frequencies in pnictides do not exceed 500 
cm -1 ,— one would be tempted to interpret the lack of 
saturation in t _1 (cj) as due to strong scattering by some 
other collective modes^ such as spin fluctuations, that 
are also the best candidates for the superconducting pair- 
ing. 

To investigate the possible role of spin fluctuations 
on the optical spectra of LaFePO we calculated the 
intraband optical conductivity of the fully interacting 
multiband model for LaFePO already introduced in 
Refs. [16.17]. In particular, the charge carriers inter- 



act with spin fluctuations modeled in terms of a Millis's 
spectrum^ 



B s{ (u>) = A a 



Tt{ujl + UJ 2 



(3) 



which emphasizes the survival of bosonic excitations up 
to large energies. Here ujq is the characteristic energy 
scale of spin fluctuations and A av is an average of the di- 
mensionless coupling constant in each of the four bands 
included in the microscopic model^i This has been 
estimated to be A av ~ 0.5 at T c = 7 K ,— to account 
for the superconducting, 21 magnetic^ and specific heat 
measurementSi 23 ' 24 This coupling is expected to be re- 
duced at T — 300 K by roughly a factor 3, A av ~ 0.17, due 
to the weakening of the antiferromagnetic correlations, 
accompanied also by a corresponding shift at higher ener- 
gies of the characteristic energy scale u>q = 60 meV ~ 480 
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The complex opti- 



cal conductivity is then computed using a finite quasi- 
particle scattering rate To = 10 meV « 81 cm -1 to ac- 
count for impurities, and an EDM is applied to these 
data to extract t~ 1 {ui), using the experimental plasma 
frequency uip = 14900 cm" M The results for ai(co) and 
t -1 (w) obtained in this way are shown in Fig. Q] as solid 
black lines, in comparison with the experimental data of 
Ref. H As one can see, while the low-energy part of 
<7i(w) is well described by the calculated intraband opti- 
cal conductivity for A av = 0.17, at ui > 500 cm -1 the the- 
oretical cti(w) deviates from the experimental data. This 
behavior is reflected also in the corresponding theoretical 
r _1 (w), which increases approximately linearly as 2A av u;, 
i.e with a slope about four times smaller than the exper- 
imental one. Indeed, within this framework, we are able 
to reproduce the experimental linear slope of r _1 (a;) at 
ui ~ 1500 cm -1 only by assuming a total coupling of the 
order of A av ~ 0.55 (red dashed line in Fig. [TJ. However 
such a large value of A av , besides being at odds with the 
dHvA and other thermodynamical measurements ) 10 ' 16 ! 7 
would result in a r _1 (w = 0), at T = 300 K, much higher 
than the experimental one. Furthermore, ai{u>), turns 
out to be too broad in the low-frequency range with re- 
spect to experimental data. The present results show 
thus that in LaFePO it is not possible to describe the op- 
tical data in the range u < 2000 cm -1 in terms of pure 
intraband transitions, even in the presence of scattering 
mediated by a spin-fluctuation spectrum with a long tail 
at high energies as in Eq. ([3]) , and regardless the specific 
value of the coupling. 

Since a model with intraband transitions alone cannot 
reproduce the optical data in pnictides, we investigate 
the possibility that the EDM analysis of the experimen- 
tal data in LaFePO is affected by interband transitions 
having a sizable spectral weight at energies comparable 
to the characteristic boson frequencies (60 meV ~ 480 
cnr 1 in our case). To elucidate in general the role of in- 
terband transitions we illustrate in Fig. [5] a few paradig- 
matic examples, where we add to the intraband term with 
A av = 0.17 from Fig. [T] (shown here as dashed black line) 
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FIG. 2: (color online) Scattering rate (panel a) resulting from 
the EDM analysis of the optical conductivity a{ui) from Eq. 
0, an d its real part o\{ui) (panel b). The intraband term 
is computed here with A av = 0.17 and the same microscopic 
parameters as in Fig. [T] The Lorentzian parameters are: (1) 
lul = 1500 cm" 1 , 7 L = 400 cm" 1 , S L = 2400 cm" 1 ; (2) 
ul = 3000 cm -1 , 7 L = 400 cm" 1 , S L = 2400 cm" 1 ; (3) 
u L = 3000 cm" 1 , j L = 8000 cm" 1 , S L = 2400 cm" 1 ; (4) 
lol = 3000 cm -1 , 7l = 8000 cm -1 , S L = 24000 cm -1 . 

a Lorentzian contribution, with a characteristic energy 
ojl, with weight S\ and width 7, which mimics inter- 
band transitions: 

a{u) = a intra (w) + f& J-. (4) 

47r w7l + ^(^i — > ) 

As one can see in Fig. [S^i, the effects on t" 1 (w) depend 
crucially on the weight of the Lorentzian and on its rela- 
tive width 7l/wl. For jl/^l <C 1 (cases 1, 2 in Fig. [2]) 
the Lorentzian peak can be clearly identified in the real 
part of the conductivity, and it also gives a sharp signa- 
ture in t -1 (u;). However, when jl/^l — 1 (cases 3, 4) 
the interband contribution adds smoothly to the intra- 
band part, see Fig. &p. This leads in turn to an overall 
smooth increase of the inverse lifetime, that is quantita- 
tively more pronounced when S"£ becomes comparable to 
the spectral weight uip of the intraband part (up = 14900 
cm" 1 as in Fig. [1]). As a consequence, in this case the 
EDM analysis of the overall conductivity would lead to 
a misinterpretation of the effects of the interaction, since 
the dramatic increase of t (uj) due to the interband 
transitions could be interpreted as due to interactions. 
As we shall argue below, this could be indeed the case in 
LaFePO. 

To establish the effective role of interband transitions 
in LaFePO we first estimate their strength and posi- 
tion by computing its complex optical conductivity using 
Density Functional Theory (DFT) in the generalized gra- 
dient approximation (GGA)?2£ We employed the full po- 
tential linearized augmented planewave (LAPW) method 
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FIG. 3: (color online) (a) Real part of the interband optical 
conductivity of LaFePO evaluated in DFT calculations, (b)- 
(d) Low-energy band-structure of LaFePO, where the thick- 
ness of the symbols represents the character jm| = 0, 1, 2 of 
each band. 



as implemented in the Wien2k codej 26 > ' 29 The Brillouin 
zone sampling for the optical calculations was done us- 
ing a very dense grid of more than 5000 k-points in the 
symmetry irreducible wedge. The complex interband op- 
tical conductivity <r(w) was evaluated using a broadening 
77 = 1 meV. We considered experimental values for both 
the lattice constants (a = 3.964 A, c = 8.512 A) and 
the internal coordinate zp = 0.3661. The interband con- 
tribution to <ti(w) is shown in panel (a) of Fig. [3J Our 
results are in excellent agreement at all energies with 
the DFT spectra in literature.— For ui < 6000 cm" 1 , the 
spectrum shows two main features: a narrow peak at 
very low energies (u « 200 — 600 cm" 1 ), and a broad 
structure starting at u) » 2500 cm" 1 . Note en passant 
that, according to the above discussion regarding the ef- 
fects of a Lorentzian interband term illustrated in Fig. [5J 
the broad high-energy feature in Fig. [3^ is expected to 
have a considerable effect on the inverse lifetime, while 
the low-energy peak is expected to produce only modest 
effects. 

The microscopic origin of the two features in the DFT 
optical conductivity is shown in the bottom panels of Fig. 
121 where we decorate the bands of LaFePO with partial 
orbital characters. In an energy range ~ 1 eV around the 
Fermi level, the bands have mainly Fe d character^ op- 
tical transitions are possible for Am = ±1, i.e. between 
to = (d z 2) and |m| = 1 (d XZl d yz ), or between |m| = 1 
and |to| = 2 (d x 2_ y 2, d xy ).- 9 The narrow low-energy peak 
is mostly due to the transitions around the F-X and Z-R 
between the three-dimensional d z i band with to = and 
the two-dimensional hole bands with character d xz , d yz 
(|m| = 1). Instead, the higher energy feature mainly 
involves transitions along T-Z between |m| = 1 states 
above the Fermi level and |m| = 2 states below it. While 
the position and shape of the first peak is very sensitive 
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TABLE I: Interband parameters of the model © expressed 
in units of cm" 1 (Sj,u^,7j) and in dimensionless units (eugh)- 



to computational details and changes for different ma- 
terials, the second, broad feature, which extends up to 
16000 cm -1 , is a robust property of the optical spectra 
of Fe pnictides and chalchogenidesi 3 ^— 

It should be remarked however that, although DFT 
calculations provide an important insight into the exis- 
tence and nature of low-energy interband transitions we 
cannot use directly the DFT data to fit the spectrum of 
Ref. 1. In fact, due to the underestimation of electronic 
correlations, DFT is known to locate the Fe d — d transi- 
tions at larger energy as compared to experiments^. For 
instance, in Ref. ■''>•''>. a direct comparison of experimental 
and DFT optical spectra of the Ba 122 compound has 
shown that the experimental position of the Fe d transi- 
tions is rescaled by a factor of two with respect to DFT 
calculations. 

Since the exact scaling factor for LaFePO is not known, 
in order to assess quantitatively the role of interband 
transitions on the spectra of Ref. |4j, we model the infrared 
interband conductivity with two Lorentzian peaks at low 
energies, 



Cr(w) = CTintraM + 



E 

i=l,2 



47T IQ^ji + i{ujf — LJ 2 ) ' 



(3) 



plus a high-energy dielectric constant ehigh to account 
for processes at higher energies. The intraband part is 
evaluated for the same microscopic parameters as in Fig. 
[1] and for A av = 0.17, which reproduces the correct ex- 
perimental behavior of the low-energy Drude-like part of 
<j\(uj) for bj < 500 cm -1 (see Fig. [I]). The Lorentzian 
parameters are meant instead to model the two groups 
of transitions shown by DFT calculations. The specific 
values for the frequency (uii), width (7^) and weight (Si) 
of the two Lorentzian contributions, chosen to fit the ex- 
perimental data for ai(uj) in the range oj G [500 : 2000] 
cm -1 , are reported in Table HI 

The result for the total optical conductivity is shown 
in black line in Fig. 3^,, where one can see the very good 
agreement with the experimental data from Ref. |4J in the 
whole range uj 6 [200 : 3000] cm -1 . The effect of the 
interband transitions on the EDM analysis is elucidated 
in Fig. UJd where we compare r _1 (o;) calculated by us- 
ing both the total optical response (dashed line) and the 
interacting intraband part only (solid line). As one can 
see in Fig. [4J3, the presence of interband transitions re- 
markably affects the magnitude and the frequency depen- 
dence of the effective EDM functions of LaFePO at very 
low energies (cj > 500 cm -1 ), in analogy with the case 
(4) of Fig. [2] discussed above. In particular, the inverse 
lifetime t~ 1 {uj) acquires a frequency dependence that is 
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FIG. 4: (color online) (a) Real part of the optical conduc- 
tivity ai(u)) computed using Eq. ([5} (red dashed line, la- 
beled as "total" in the caption), along with the intraband 
part (solid black line) and the interband part (solid green 
line). The symbols represent the experimental data of Ref. 
4. (b) r _1 (o;) extracted from the extended-Drude-model us- 
ing the total conductivity (dashed red line) or the intraband 
part only (solid black line), along with the experimental data 
(symbols) (c) Approximate electron-boson Eliashberg func- 
tion W(u)) extracted from t~ 1 (ui) with the second-derivative 
method (see text) using either the total conductivity (dashed 
red line) or the intraband part only (solid black line). 



much stronger than the one of the intraband contribution 
alone, with a linear slope which is w 4 times larger than 
the one arising from the actual intraband contribution. 

In the spirit of the EDM, the frequency dependence 
of r _1 (o;) is often used to extract the underlying spec- 
trum a 2 F(ui) of the retarded interaction supposed to 
generate it. A crude approximation to obtain the in- 
teraction spectrum is the second derivative method, in 
which a 2 F(uj) is approximated by W(ui), where W(uj) = 
(l/2Tr)d 2 [ujT- 1 (w) / dw 2 ]. 35 ' 36 This method is usually suf- 
ficient to capture the order of magnitude of the spectrum, 
although not its fine details. The W(u>) corresponding to 
the total and to the pure intraband T _1 (tj) are shown 
as dashed red and solid black lines, respectively, in Fig. 
Hfc, showing that the spectrum extracted from the total 
r _1 (w) is significantly larger than the other one. In- 
deed, the "effective" coupling A g = 2 f^" dwW (uj) / uj , 
with u) c = 3000 cm -1 , obtained from the intraband term 
results A c ff = 0.27, larger but of the same order of mag- 
nitude as the microscopic value used here, A av = 0.17, 37 
whereas the effective coupling estimated using the total 
t -1 (oj) would give A c ff = 1.31, significantly larger than 
A av . This observation can thus explain why a direct EDM 
analysis of optical data, without a proper subtraction of 
the interband processes, gives much larger values of the 
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coupling^ than what estimated by other probes i 16 i 23 i 24 
In summary, we have shown that the EDM analysis in 
pnictide systems can be strongly biased by the presence 
of significant interband contributions in the mid-infrared 
region of the optical spectrum. In particular, the strong 
increase of the inverse lifetime r _1 (w) , observed in several 
compounds, cannot be straightforwardly attributed to 
strong-coupling effects, as it would be the case for opti- 
cal spectra dominated by intraband interactions. Even 
though we focused explicitly on LaFePO, we expect sim- 
ilar effects to hold in other classes of pnictides as well, 
due to the very similar band structure. In particular, 
the presence of low-energy interband transitions can pro- 
vide an alternative explanation for the optical spectra of 
122 compounds, where the flat region of the <7i(w) has 
been usually fitted with a Drude component having an 
unrealistically large scattering rate.— ^ 
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